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Cyclotron resonance in a two-dimensional electron gas with long-range randomness 
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We show that the the cyclotron resonance in a two-dimensional electron gas has non-trivial prop- 
erties if the correlation length of the disorder is larger than the de Broglie wavelength: (a) the 
lineshape assumes three different forms in strong, intermediate, and weak magnetic fields (b) at 
the transition from the intermediate to the weak fields the linewidth suddenly collapses due to an 
explosive growth in the fraction of electrons with a diffusive-type dynamics. 



-a 
c 
o 

(N 
> 

in 



o 

oo 
On 



i 

C 

o 
o 



X 



The cyclotron resonance (CR) is one of the basic tools 
for studying the electronic properties of physical systems 
in an external magnetic field. A very interesting exam- 
ple of such a system is a two-dimensional electron gas 
(2DEG), which exhibits the quantized Hall effect if the 
magnetic field is sufficiently strong. The CR can be stud- 
ied by measuring the transmission of an electromagnetic 
signal of some frequency u> through the 2DEG. If the 
absorption is weak, the change in the transmission is 
proportional to the real part of the dynamical conduc- 
tivity cr xx (uj), which is the average over the active and 
inactive circular polarizations, Kea xx (uj) = [Recr+(oj) + 
Recr_(aj)]/2 (for linearly polarized signals). The active 
polarization's contribution Reo+(o;) has a peak at u) 
close to the cyclotron frequency ui c of the external mag- 
netic field. The disorder-related zero-temperature width 
and shape of this peak are the subjects of this Letter. 

The theoretical study of the CR in the 2DEG has been 
initiated by Ando |j| . Although many pieces of the com- 
plete picture has been put in place by him and later by 
other authors |2|J^| , the consistent description of this phe- 
nomenon exists only for the case where the correlation 
length d of the random potential U(x,y) acting on the 
electrons in the disordered 2DEG is smaller than the de 
Broglie wavelength kp 1 (k F is the Fermi wavevector in 
zero magnetic field). In this case the effect of the random 
potential is described by a single quantity, the transport 
(or momentum relaxation) time r. It was understood, 
however, that the conventional Drudc-Lorcntz formula 
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(<7o is the zero- field conductivity), which gives the 
Lorentzian peak with the half width at half maximum 
(HWHM) of t -1 , applies only in the uninteresting case 
uj c t <C 1. In the other limit (uj c t 3> 1) the CR lineshape 
is non-Lorentzian and has a much larger width Q| 



due to the formation of discrete Landau levels. This be- 
havior of Awx/a is illustrated by the thin line in Fig. [l| 
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Awi/2 = 0.73 y/ujc/i 
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FIG. 1. Dependence of the CR linewidth on (oj c t) -1 , 
the quantity inversely proportional to the magnetic field. 
Thick line: our results for the long-range potential with 
given W/E F and r [for the case k F d > (E F /W) 2/3 ]. La- 
bels (0), ((?]), and (g) correspond to the equation num- 
bers. Thin line: short-range potential with the same r. 
Dash-dotted line: the width v of an additional structure 
near the very resonance (see text). 

Here and below Awi/2 is the median width defined by 

<-Ao>i/ 2 y r°° 

/ divRecr+(uj + lo c ) = — j dui Re a+(uj + uj c ). 
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Using the median width instead of the conventional 
HWHM is more adequate because the CR lineshape can 
be rather intricate for a long-range random potential. 

For simplicity, we will consider a model of a Gaussian 
random potential whose correlator decays sufficiently fast 
at r > d and does not possess any other characteristic 
scales besides d. The rms amplitude of U will be de- 
noted by W. Until the very end we will assume that 
k F d 3> (E F /W) 2 / 3 , where E F is the Fermi energy. As 
one can see from Fig. [l] illustrating our results for this 
case, the dependence of Aa>i/2 on the magnetic field is 
non-monotonic. Even more remarkable, Ao;x/2 exhibits 
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a rapid collapse to its classical value of r 1 in the vicin- 
ity of the point lo c t ~ (E F /W) 2/3 > 1. The derivation 
of these results is based on the picture of the "classical 
localization" which we will now discuss. 

If the random potential U is smooth and the time scale 
on which we study the motion of an electron is not too 
long, then it can be described classically, as a motion 
of a single particle with energy E F . (We neglect the 
interaction). If the magnetic field is not too low, the 
motion can be decomposed into a fast cyclotron gyra- 
tion and a slow drift of the guiding center p = (p Xl p y ) 
of the cyclotron orbit. For magnetic field in the nega- 
tive ^-direction, p x — x + (v y /ui c ) and p y = y — (v x /ui c ), 
where x and y are the coordinates of the electron and 
v = — vp (sin.0, cos9) is its velocity. We will call the 
magnetic field strong if the cyclotron radius R c — vf/u> c 
is smaller than d. It is usually assumed that in such 
strong fields the guiding center drifts along a level line 
U(p) = const of the random potential, which is typi- 
cally a closed loop of size d. Recently it has been re- 
alized [QJ|] that the drift approximation is also valid in 
the intermediate field regime 1 < R c /d < (E F /W) 2 ^ 3 
[the same as (E F /W) 4 ^ 3 < uj c t < (E F /W) 2 because 
t ~ (d/vp)(Ep/W) 2 ]. The point is that in this regime 
the cyclotron gyration is still sufficiently fast, so that the 
guiding center remains practically "frozen" during one 
cyclotron period. Therefore, the guiding center motion 
is determined by Uq, the random potential averaged over 
the cyclotron orbit, 

f 2rT dO 

U (p,R c ) = / —U(p x + R c cos9,p y + R c sm6). 
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The guiding center is still bound to one of the level lines, 
but those are the level lines of Uq, not U jf|. The rms 
amplitude Wo of the potential Uo is by a factor of y/d/R c 
smaller than W due to the averaging, but the correlation 
length is the same; therefore, a typical level line of Uo 
is still a loop of size d. The frequency uid of the guiding 
center motion along such a loop (the drift frequency) is 
given by Wq/ (muj c d 2 ), which can be cast into the form 



u d = (R c /d)y / u c /T, R c > d. 



(3) 



Obviously, the electrons on the periodic orbits are (clas- 
sically) localized and do not participate in the DC trans- 
port. However, a more accurate analysis || reveals that 
a very small fraction of the order of e~^ c / Ud ^ 2/3 of the 
trajectories remains delocalized. Such trajectories form 
a stochastic web in the vicinity of the percolation con- 
tour. The stochastic web rapidly grows with decreasing 
magnetic field and turns into a stochastic sea at 



R c /d = {E F /Wf/z, 



(4) 



where Ud — lo c - In even lower magnetic fields (the weak 
field regime) the stochastic sea spans almost the entire 



phase space. Correspondingly, the static conductivity is 
exponentially small, <x xx (0) cx e"^"^^ /3 when u>d <C uj c 
(the strong and the intermediate field regimes), rapidly 
blows up near u>d — ui c point (the boundary of the inter- 
mediate and the weak field regimes), and finally crosses 
over to the Drude-Lorentz formula ([l]) in weak fields || . 

Our goal is to show that the dynamical conductivity 
also exhibits a rapid change near the u>d — lo c point: the 
aforementioned collapse of the CR linewidth. This can 
be done on the basis of the classical formula for a+(u>), 
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Suppose that Alu = ui — u) c is smaller than ui c by abso- 
lute value, then the full equation of motion for 9 can be 
replaced by its average over the cyclotron period, 



dt 



Suj, Slu 
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d 



mui c R dR 



U (p,R) 



(6) 



R=Rc 



where Suj is the local correction to the gyration frequency 
of the velocity vector. The origin of such a correction is 
quite clear. Indeed, the random potential creates an ad- 
ditional centripetal force, which is either parallel to the 
Lorentz force and thus speeds up the cyclotron gyration, 
or anti-parallel, which slows it down. 

In the strong and intermediate field regimes most of 
the guiding centers are localized on the periodic orbits of 
small size < d. The correction 5ui does not vary much on 
this length scale. In fact, to find Awi/2 up to a numerical 
factor, we can use the approximation Su> = const. Sub- 
stituting this into Eqs. (||) and @, we find that Recr + (aj) 
has the form of a Gaussian peak of width [&j] rms , and so 



W 

Awi/2 ~ [cMrms = const -r , R c < d, (7) 
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(8) 



VWH, Rc > d. 



The broadening of the CR line is of the inhomogeneous 
type, and comes from the places where Aui = Suj |Q. 
Equations (0) and (^) show that linearly in- 

creases as a function of (uj c t)~ 1 in the strong-field regime, 
reaches its maximum at R c ~ d, and then decreases in 
the intermediate fields, see Fig. [j]. 

The situation is different in the weak fields where most 
of the trajectories are extended and ergodic. In this case 
the broadening is of the homogeneous type and is much 
smaller than [(5w] rms because of the motional narrowing. 
The crossover between the two types of broadening oc- 
curs at the low-field end of the intermediate-field regime, 
which we study below. (However, a broader physical dis- 
cussion will be resumed in the concluding remarks). 

According to Eqs. (||) and (§h, [<W] rms <C oJd in this 
regime. Because of this inequality, the orbits of low 
frequency Q <C UJd whose perimeter length is typically 



2 



(u>g/fl.)d turn out to be important. To show this we 
present o + (w) as follows, 



ct + (w) = 



<T0 



dnf(n)i(n,Au), 



I(CL, Aw) = J dt (exp iAujt - i J 



iAcot-i I dt'Suj(t') 
o 



(9) 
(10) 



where (• ■ -)n denotes the average over the orbits of fre- 
quency f2, and f(£l) is the probability density of finding 
an electron with energy Ep on such an orbit (with the 
convention £1 = +0 for the unbound trajectories of the 
stochastic web). The correlations in Sco(p) decay suffi- 
ciently fast (~ p -3 ) with distance; hence, the high-order 
correlators of 5oj can be neglected, and we arrive at 
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dt exp 



iAtot - 



dt 



dt 2 (6u)(ti)6u)(t 2 ))n 



For the same reason we can approximate the second-order 
correlator by the sum of isolated short pulses, 



the role of the effective transport time, and so Eq. jl5| ) 
is expected by analogy with Eq. ([!]) . 

To finish the calculation of a+(w) we need to know 
/(f2). It is determined by the statistical properties of the 
potential U, viz., by the exponent H in the equation 



2q 



-2H 



dqq.Jg(qR c )C(q) cx q 



where Jo is the Bessel function and C(q) is the Fourier 
transform of C(r) = (U(0)U(r)). In a practically rele- 
vant case where the potential is created by randomly po- 
sitioned donors confined to a narrow layer separated by 
the distance d from the 2DEG, C{q) = C(0)e- 2qd , so that 
H = -1 if q < R- 1 and H = -1/2 if R- 1 < q < cT 1 . 
In both cases 

f(fi) ~ sn- 1 (n/u, d y, n wch <n<u d , (i6a) 

~ (Q wcb /uj d ) s S(Q - 0) at other fi, (16b) 



(M*0M*a))n= E P(h-t 2 -^), (11) 
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9. 



In 



where P(0) = {5lj 2 ) and P(t) < P(0) for \t\ > u 
fact, if both f2 and |Aw| are much smaller than uid, the 
actual functional form of the pulses is unimportant, and 
we can approximate P(t) by 2v5(t), where 



v ~ (Suj 2 )/uj d ~ w c /(w d r). 



(12) 



After such approximations and some algebra we find: 



Re J 



^ c n exp 



n— — oo 



c n = I dx cos(wnx) exp 
Jo 



tt(Aw - nVi) 2 

itv{2x — x 2 ) 
2Q 



(13) 
(14) 



If 3> v, then Co ~ 1 and c n ~ i//(irn 2 fl), n ^ 0, so 
that Re/(f2, Aw) as a function of Aw is the sum of nar- 
row Gaussians centered at equidistant points Aw = nQ. 
In the opposite limit the Gaussians merge into a single 
Lorentzian, 



Re 1(0, Aw) = — 
v (A 



(Aw) 2 + . 



f2 < v, |Aw| < w d . (15) 



There is a simple way to understand this equation. As 
discussed above, the random potential modifies the lo- 
cal gyration frequency of the velocity vector by Slo. The 
sign of the correction changes randomly each time the 
guiding center travels the distance ~ d. This causes the 
angular diffusion of the velocity vector in addition to the 
regular precession with frequency w c and therefore, the 
momentum relaxation. The time uj d /(Suj 2 ) = v^ 1 plays 



where s is very close to 1/2 [see Isichenko || for de- 
tails; f2/(f2) corresponds to his F(Q,)]. The quantity 
f2 wo b ~ LOde^^"/^^ 2 3 is the drift frequency of such or- 
bits within the stochastic web that give the dominant 
contribution [Q to the DC transport. 

We can now substitute Eqs. ([0§{L§) into Eq. (|). bmce 
s < 1, only the n — term in Eq. (|l^) is impor- 
tant for v < fl < u! d - Also, it is safe to assume that 
u d /v < e 2 /^ 8 " 1 ) §. With this in mind, we obtain 
Re<7 + (w) = (ctq/vt) [Si (w) + S^w)], where 




In 



w d i + {Abj/vy 



max {(Aw) 2 , tt c i>} 



, |Aw| < y/uJdV, 



fl c = max{f2 we b, v}. It is easy to see that as long as 
^web < v <C w<j, Awi/2 is determined by ^(w), and thus 
is of the order of ^JuJdP ~ v/ w c /t, in accordance with 
Eq. (@). The CR lineshape for this case is depicted in 
Fig. ||a. 
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FIG. 2. The CR lineshape (schematically), (a) before 
the collapse (b) in the course of the collapse (c) immedi- 
ately after the collapse. 
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Let us now discuss the aforementioned collapse of the 
CR linewidth. il wc b grows exponentially as the mag- 
netic field decreases and eventually becomes larger than 
v. From this moment, the maximum value of Si rapidly 
decreases and that of S\ does the opposite. It gives rise 
to a narrow Lorentzian peak at the center of the CR line 
(see Fig. |^b). At the tod = uj c point where f2 wc b ~ 0Jd, it 
is this peak that determines Au>x/2- Thus, Au>i/2 drops 
to the value of v = uj c /(ujdT) — r _1 (Fig. ||c). From this 
point on, i.e., in the weak- field regime, Eq. ([!]) applies. 

Several comments are in order here. Formula (Q) for 
the strong-field regime is not new (see Refs. |lj-[|). How- 
ever, its derivation has remained unsatisfactory: Ando's 
work |jj is based on the self-consistent Born approxima- 
tion (SCBA), which is invalid in this regime; Prasad and 
Fujita ignored the localized nature of the particle mo- 
tion; Bychkov and Iordanskii || de facto assumed that 
the electron's trajectory (in the absence of interactions) 
is Sui(p) = const instead of the correct Uo(p) — const. 

So far we have studied the case kpd ^> (Ep/W) 2 ^ 3 . In 
the opposite limit an additional "quantum" region ap- 
pears on the phase diagram of Fig. | p0|| . Its bou ndary 
is formed by the lines I = d where I = ^h/muj c is the 
magnetic length, and w c r 9 = 1 where r q ~ t /(kpd) 2 is 
the quantum lifetime. The lower edge of this region, i.e., 
the line kpd = 1 was discussed in the beginning of this 
Letter. Now we will show that the collapse of the CR 
linewidth occurs in the case 1 -C kpd <C (Ep/W) 2 ^ 3 as 
well, and that its position on the phase diagram is given 
by the line I = d (which is the same as R c /d = kpd). 
Contrary to the "classical" case, the quantum collapse 
takes place when the Landau levels are well separated. 




FIG. 3. Phase diagram of the CR (double log-scale). 
iV is the index of the Landau level closest to Ef- The 
bold line marks the location of the linewidth collapse. The 
paths A and B, which the 2DEG follows as the magnetic 
field decreases, correspond to the classical (Fig. [l]) and the 
quantum (Fig. ^) cases, respectively. 

The CR lineshape in the quantum region can be stud- 
ied within the SCBA. Doing the integration in Eq. (2.12) 
of Ref . jt) for the case of Ep positioned between the Lan- 
dau levels and |Aw| <C y/us c /T q , we find that the CR line 
is roughly Lorentzian, while its median width is given by 

Au> 1/2 = OMr' 1 ^^. (17) 



On the other hand, in the intermediate- field region Awi/2 
is given by Eq. (Eh. Therefore, the crossing of the I = d 
line causes the drop of Au>i/ 2 by a large factor of kpd 
(Fig. ||) . The physics of this "quantum" collapse is quite 
similar to that of the classical one: the crossing of the 
line I = d is accompanied by an explosive growth of the 
quantum localization length of the states near the Lan- 
dau level center (see Jl0| for discussion) . The states with 
sufficiently large localization length contribute to o~+(uj) 
in the form similar to Eq. (|l|) but with r replaced by 
the effective transport time. Since the density of states 
at the Landau level center is larger than the zero-field 
density of states by a factor of ^/u> c T q Q , this effective 
transport time is smaller than r by the same factor in 
agreement with Eq. ([l7|). 




kfdiW/Epf (kpdf 1 (lo c t) A 
FIG. 4. Same as Fig. | but for the path B of Fig. |. 

Until now we considered W a fixed parameter of the 
theory. In a more realistic model (see above) W is 
determined by the concentration m of randomly posi- 
tioned donors and the screening properties of the 2DEG. 
Away from the strong-field regime, the screening is 
very much the same as in zero field, and one obtains 
W ~ Epy/ni/(k F d). Thus, the "quantum" case kpd <C 
(Ep /W) 2 / 3 is realized if rii is sufficiently low, <C kp/d. 
On the other hand, in the strong-field regime (far from 
the collapse), the amount of screening, W, and Awi/ 2 os- 
cillate with the filling factor |llj . A consistent description 
of such oscillations must take into account the nonlinear 
screening effects ]l2] ]. 

In conclusion, let us discuss the experimental side. The 
non-monotonic dependence of Awx/2 on t ne magnetic 
field with the maximum at R c ~ d has been observed 
in Ref. [ fL3f , in agreement with our theory. Arguably, the 
collapse of the CR linewidth has also been seen (for the 
400A spacer sample). However, a decisive confirmation 
of the latter prediction requires further experiments. 
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